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H , We investigate the cosmic variance of the skewness of the cosmic microwave back- 

ground (CMB) anisotropies in an inflationary model which leads to the baryon isocur- 
vature scenario for the cosmic structure formation. In this model, the baryon number 
fluctuations are given by a sinusoidal function of a random Gaussian field. We find that 
the skewness is very small in comparison with that of the fluctuations which obey Gaussian 
statistics. 



Abstract 



The existence of the nearly isotropic cosmic microwave background (CMB) is a firm ev- 
idence of the hot big bang cosmology. At the same time, it is believed that its anisotropies 
carry valuable information about the early history of our universe. Recently, the COs- 
mic Background Explorer (COBE) satellite team reported the detection of temperature 
anisotropies in the CMB on scales greater than 10° (Smoot et al. 1992). These fluctua- 
tions are particularly important because they directly reflect the gravitational potential 
fluctuations at z ~ 1000. They are hence tightly related to the structure formation of the 
universe, and their data are powerful clues to know the origin of the cosmic structure. 

The angular two-point correlation function is commonly used in a statistical analysis 
of the temperature fluctuations. However, its knowledge alone cannot distinguish the 
statistical properties of the fluctuations. To do so, one needs to know at least the next 
order correlation, i,e, the angular three-point correlation. In the case of the temperature 
fluctuations which originate from the adiabatic curvature perturbations in the inflationary 
cosmology, it is shown that its imprint on the angular three-point correlation function is 
very small (Falk et al. 1993; Gangui et al. 1993), and that the contribution from the 
cosmic variance will dominate over the signal (Srednicki 1993; Scaramella 1991). 

Recently the COBE team reported the results of the data analysis of the three-point 
correlation function based on the DMR flrst-year sky map (Hinshaw et al. 1993). In 
terms of a Monte Carlo analysis with an assumed level of the instrumental noise, they 
found an evidence of the non- vanishing three-point correlation function in the data which 
is consistent with the prediction of the standard adiabatic cold dark matter (CDM) model 
based on the inflationary cosmology. However, it is also true that the low signal-to-noise 
level in the data prevented them to reproduce the three-point correlation function in a 
clear form, and hence they could only place an upper limit on the amplitude of it as a 
flrm conclusion. According to them, further several-year observations will produce clear 
data. 

In this situation, it is worthwhile to consider the angular three-point correlation in 
other possible scenarios of the structure formation besides the standard CDM scenario. 
Among them, a viable alternative is the Peebles isocurvature baryon (PIB) scenario of the 
structure formation, in which the universe consists of only baryonic matter and radiation 
but no non-baryonic dark matter (Peebles 1987), based on the observation Qq ~ 0.1 
(Peebles 1986; Tyson 1988). 

Furthermore, the second year COBE-DMR data reported recently indicate a power-law 
index of the primordial density fluctuations which is higher than that of the Harrison- 
Zeldovich spectrum (Bennett et al. 1994). In addition, several other recent CMB ex- 
periments on degree scales have reported the detection of a high amplitude temperature 
anisotropy (Gundersen et al. 1993; Hancock et al. 1994), which may be regarded as an 
evidence in favor of the PIB scenario (Hu & Sugiyama 1994), though not all of the results 
of degree scale experiments seem to be consistent, which may be due to non-Gaussian 
statistics (Luo 1994). 

However, the shortcoming of the PIB scenario, in the theoretical sense, is that it 
assumes a very ad hoc spectrum of the primordial density fluctuations to explain the 
structure formation and to satisfy the observed isotropy of the CMB on large angular 
scales, namely isocurvature baryon fluctuations with a steep spectral index. 

Now, an interesting microscopic mechanism to produce the seemingly ad hoc spectrum 



for the PIB scenario has been proposed (Sasaki & Yokoyama 1991; Yokoyama & Suto 1991) 
In this mechanism, the soft CP violation is induced by a spatially-varying Majoron field, 
A{x.), associated with a heavy Majorana lepton field which decays three quarks or three 
anti-quarks, and the space-dependent net baryon number, 

i?(x) = B, sin (^) (1) 

is generated, where / is the mass scale of the Majoron field and -B* is a constant determined 
by the coupling constants of the model. In particular, Sasaki and Yokoyama (1991) found 
that a power-law inflation model provides the most appropriate class of power spectra for 
the PIB scenario. In this case there appears a characteristic scale kc, which is determined 
by the model of inflation and the mass scale /, and the power spectrum approaches the 
white noise as k/kc —>■ and almost the scale invariant one as k/kc — > oo. Thus it has a 
nice feature that the amplitude of small scale fluctuations does not become too large on 
very small scales. Further, they have shown that there exists a natural particle physics 
model that can provide an appropriate initial condition for the PIB scenario, implemented 
in power-law inflation with the power-law index around n^ ^ 10 ~ 20. The only difference 
from the original PIB scenario is that one naturally expect the presence of a cosmological 
constant that makes the universe spatially flat in this inflationary isocurvature baryon 
model. Detailed statistical properties of the baryon fluctuations in this model has been 
studied by Yamamoto et al. (1992). Considering the viability of the PIB scenario, it is 
important to examine if this model for the PIB scenario has statistical properties which 
are observationally testable. 

In this letter, based on the results of Yamamoto et al. (1992), we consider the skew- 
ness of the CMB temperature fluctuations in this scenario. We denote the temperature 
fluctuation field by ^(x, 71), where x specifies the position of the observer, the unit 
vector 7i points a given direction from x. The temperature fluctuation can be evaluated 
using the isocurvature perturbation theory (Kodama & Sasaki 1986). The baryon number 
fluctuations correspond to entropy perturbations which give rise to gravitational poten- 
tial perturbations after the universe becomes matter-dominated. Then the temperature 
fluctuations directly trace the primordial baryon number fluctuations on large angular 
scales, 

^(x,7i)= ci?(xi), (2) 

where c is a proportional constant (Sasaki & Yokoyama 1991), Xi := x + ro7i, Vq = 

H'^ Jo dy/^^oy + (1 - ^o)y^ ~ 2VLq^-^Hq^ (^o ~ O.l), and Hq is the Hubble constant. 
Thus the statistics of the temperature fluctuation can be understood through that of 
S(x). 

Now we focus on the skewness of the CMB anisotropy. 
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where dVt^^ denotes the integration over the direction 71. As is clear from Eqs.(|^) and (^, 
we have (5*) = 0, where ( ) denotes the ensemble average on the position of observer x. 



Then we consider the cosmic variance of the skewness, which can be written as 
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The ensemble average of the temperature fluctuations can be replaced by the correlation 
of the baryon number fluctuations by using Eq.(|^). The random field -B(x) follows a 
peculiar statistic, because it is a sinusoidal function of the Gaussian random field v4(x). 
However, we can evaluate the 2m-point correlation function of -B(x) by using the following 
formula that is derived in by Yamamoto et al. (1992), 
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where (3 := H^/Sn f , Hj is the Hubble parameter at the time when the inflation ends, Up 
is the power-law index of the power-law inflation, o"j (z = 1, 2, ■ ■ ■ , 2m) takes the value ibl, 
and Yl'a means that the summation is taken all over the combinations (o"i, cr2, ■ ■ ■ , (J2m) 
with J2'j=i o^j = 0- After straightforward but tedious calculations, we find. 
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Note the appearance of the 9th power of the two-point correlation in the first term. This 
term arises because there are 3x3 independent combinations between Xj (z = 1,2,3) and 
Xj (j = 4,5,6) before we take the coincidence limits. The two-point correlation of -B(x) 
is given by (Sasaki & Yokoyama 1991) 
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where 
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of the multipole moments as 



The ensemble average of the angular two-point correlation function is given by C{a) : = 
^(x, 7i)^(x, 72)) = c^(i?(xi)_B(x2)), where cosa = 71 ■ 72. This is written in terms 
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Here the window function Wi = exp{—l{l + l)/2cr^), a = 17.8, is inserted in order to 
compare it with the result of COBE (Smoot et al. 1992). 
Then, from Eqs.(^),(|) and (|^), we obtain 
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where we have used the fact that the monopole term (and also the dipole term) in the 
temperature fluctuations should be subtracted from the data. The variance of the skew- 
ness in the Gaussian statistics takes a similar form (Srednicki 1993), but differs from the 
present case in the powers of C{a). 

For comparison with other papers we normalize the variance of skewness by C{a = 0)^, 
and consider the normalized root mean square skewness, i.e., Sr.m.s. '■= (5'^) /C(0)^/^. 
We can show that this can be written as 

Sr.m.s. = (^^e^-^^A, (11) 

where 
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For the present scenario to be successful, we assume Up/S < 1 (Sasaki & Yokoyama 1991). 
Hence, the coefficient of A in the Eq. (pUD is of the order of unity. Therefore we focus on 
the numerical factor A. If the two parameters, i.e., np (the exponent of the power-law 
inflation) and Xc '.= r^k^ (the ratio of the present horizon size tq to the characteristic 
scale 1/fcc) are specified, we can calculate it numerically. The Table 1 shows the results 
of numerical integration of A for various values of Xc in the case Up = 10 and the Table 2 
does in the case Up = 20. 

Once Hp and Xc are fixed, we can also obtain the angular two-point correlation function 
C{a) from Eq.(^. We numerically fitted this angular two-point correlation function C{a) 
normalized by C(0) to that due to the density fiuctuation with the power-law spectrum 

((Ap/pj ) oc k'^"^ , and calculated the best fitted value of ries. Thus ries is the effective 
spectral index of the fluctuations on the COBE scale; ries ^ 2 for Xc ^ 1 while n^s ~ 1 
for Xc ~ 1. 

Let us roughly evaluate an upper limit on A. In the case Xc ~ 1, for which the 
spectrum is almost scale invariant, we have di ^ 7r/2/(/ + 1)?^^. Then we flnd A ~ O(np^). 
However, this case is not appropriate for the PIB scenario. When Xc ^ 1, which is 
the case of our interest, the value of A is very small as shown in the Tables 1 and 
2. Thus the cosmic variance of the skewness turns out to be very small in comparison 
with the one obtained for the Gaussian statistics with the Harrison-Zeldovich spectrum, 
Sr.m.s. — 0.18 (Srednicki 1993). This suggests that the cosmic variance of the angular 
three-point correlation function is also very small. We suspect that the higher order 
correlation is always suppressed in this model. 

If we accept the data analysis of the COBE team (Hinshaw et al. 1993), i.e., the 
detection of the nonvanishing three-point correlation as they claim, it is difficult to explain 
it within the present model of the PIB scenario. However, either the actual amplitude of 
the skewness or the form of the three-point correlation function has not been obtained due 



to the high noise leveL According to them, the noise level of the three-point correlation 
function diminishes in proportion to {time)~^^'^ as the data accumulate. It is necessary 
for us to wait for a few more years before a definite conclusion may be drawn. 
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Table 1 
Skewness and effective spectral index for Up = 10. 



Xe 


A 


^cff 


10^ 


3.2 X 10-13 


2.6 


10^ 


3.0 X 10"9 


2.1 


10^ 


1.4 X 10-"^ 


1.7 


103 


5.5 X 10-5 


1.4 


102 


3.6 X 10-^ 


1.17 


10 


7.3 X 10-^ 


1.04 



Table 2 
Skewness and effective spectral index for n^ = 20. 



X, 


A 


ncff 


lOii 


2.8 X 10-1^ 


1.9 


10^ 


6.5 X 10-1° 


1.6 


10^ 


3.7 X 10-^ 


1.33 


10^ 


1.0 X 10-5 


1.17 


103 


5.6 X 10-5 


1.07 


10 


9.6 X 10-5 


1.01 



